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Quantum walks, both discrete (coined) and continuous time, on a general graph of N vertices 
with undirected edges are reviewed in some detail. The resource requirements for implementing a 
quantum walk as a program on a quantum computer are compared and found to be very similar for 
both discrete and continuous time walks. The role of the oracle, and how it changes if more prior 
information about the graph is available, is also discussed. 
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I. INTRODUCTION 

One of the most important tasks on the theoretical side 
of quantum computing is the creation and understanding 
of quantum algorithms. The recent presentation of two 
quantum algorithms based on quantum versions of ran- 
dom walks is particularly important in this respect, since 
they provide a new type of algorithm which can show an 
exponential speed up over classical algorithms, to add to 
those based on the quantum Fourier transform. Childs 
et al. Q have produced a scheme for a continuous time 
quantum walk that can find its way across a particular 
graph exponentially faster than any classical algorithm, 
while Shenvi et al. proved that a discrete quantum 
walk can reproduce the quadratically faster search times 
found with Grover's algorithm for finding a marked item 
in an unsorted database. For an overview of the develop- 
ment of quantum walks for quantum computing, see the 
recent review by Kempe j3j- 

The importance of classical random walk techniques 
for computational physics motivates the study of quan- 
tum walk algorithms both for understanding in principle 
the power of quantum computation, and for the possi- 
bility of a practical quantum computer that can solve 
problems more efficiently. It is important to make a fair 
comparison between quantum and classical algorithms 
when determining whether the quantum method is re- 
ally faster. Resource counting is difficult at the best of 
times, and for the quantum walk algorithms presented so 
far, the presence of an oracle is an added complication. 
One point that has not been made explicit by most au- 
thors to date is that if, as is usually the case, the outcome 
of the quantum walk will be obtained by measuring the 
position, the quantum walk algorithm should have the 
position space encoded into a qubit register (or equiva- 
lent Q). This is simply because for a classical random 
walk algorithm, the location of the particle can be en- 
coded in a binary string, so we will penalize ourselves by 
requiring exponentially more resources if we adopt unary 
encoding for our quantum walk. Childs et al. |l| explic- 
itly perform this encoding, but without remarking on this 
point, since their primary motivation is to show that they 
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can implement a continuous time Hamiltonian evolution 
efficiently on a discrete quantum circuit. On the other 
hand, this places the proposed physical implementations 
of a quantum walk |a,|6|jl3 firmly in the realm of physics, 
as all of them have the position space of the walk set up 
in such a way that it cannot be measured as a binary 
encoded bit string. None of them claim they are trying 
to perform algorithms with these walks, and the exper- 
iments are of considerable interest in their own right as 
examples of coherent control in physical systems. For a 
discussion of the fundamental characteristics of quantum 
walks in physical systems, see Ref. 0. 

There are important open questions surrounding the 
relationship between the two main formulations of q uan- 
tum walks in the literature, discrete time |!| Efl UE [CI j 
and continuous time Classically, the discrete and 

continuous versions of random walks can be related in 
a straightforward manner by taking the limit of the dis- 
crete walk as the size of the time step goes to zero. In 
the quantum case, the way the systems are commonly for- 
mulated, the discrete time and continuous time versions 
have different sized Hilbert spaces so one cannot be the 
limit of the other without showing that the larger system 
is actually restricted to a subspace the same size as the 
smaller under suitable conditions. Also open, but gener- 
ally expected to be true, is whether discrete and contin- 
uous time quantum walks have the same computational 
power. They give broadly similar results for algorithmic 
properties such as mixing times and hitting times in all 
cases where both forms have been applied to the same 
problem though their detailed dynamics are different. 

The main aim of this paper is to describe a quantum 
walk process for both discrete (coined) and continuous 
time quantum walks on a very general graph side by side, 
so they can easily be compared and contrasted. The gen- 
eral case of the discrete time quantum walk has only ap- 
peared in the literature in Watrous 12], in the form of 
a unitary process on a directed graph derived from the 
undirected graph under study. This formulation is ex- 
actly equivalent to using a coin combined with the origi- 
nal (undirected) graph, but physicists seem to prefer the 
coined version, so it is useful to write out the coined 
quantum walk in the general case. We then observe that 
when a continuous time quantum walk is implemented 
as a quantum computer program, the structure of that 
program turns out to be very similar to the discrete time 
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walk on the same graph. 



II. CLASSICAL RANDOM WALK ON 
GENERAL GRAPH 

We start by defining the type of graph we will consider 
here. A graph G(V, E) has a set of N vertices v € V and 
a set of edges e € E joining pairs of vertices. The edges 
are undirected, that is, e xy = e yx for the edge joining 
vertices v x and v y . We assume there is at most one edge 
between any pair of vertices. The maximum number of 
edges that can meet at each vertex - the degree of the 
graph - is thus (N — 1), and the maximum number of 
edges we could have in total is N(N — l)/2, one between 
every possible pair of vertices. The degree of a vertex is 
the number of edges that meet at it, and can be zero if 
there are no edges joining that vertex to any other. 

The description of the graph is conveniently expressed 
as the adjacency matrix A, which is of dimension N x N, 
symmetric, with element A xy = 1 iff e xy G E and zero 
otherwise. Since there are N(N — l)/2 possible edges, 
A actually contains two copies of this information, but 
that duplication allows us to use a matrix rather than 
a triangular table, and avoids having to first determine 
whether x > y or vice versa when reading out values of 
A xy . The ordered way in which information is stored 
in A is an example of a hash table that allows any bit 
within it to be accessed with a single query specifying the 
values of x and y. The full information describing a single 
edge needs to be given not by a single bit, but by the set 
(x, y, A xy ) that identifies the edge by its vertex indices. 
With no other prior information about the graph, the 
adjacency matrix is the most efficient way to store and 
access the description. 

Already our general graph has several conditions at- 
tached to it, undirected edges and only one edge be- 
tween any pair of vertices. Graph theory is such a broad 
field that variants suggest themselves almost anywhere 
you look. Directed edges (that can be traversed in only 
one direction) give a completely different character to 
the graph, like the difference between reversible and ir- 
reversible processes. They are beyond the scope of this 
paper (for reasons which will become clear later), the 
reader can refer, for example, to Severini [LH llfij where 
a (different formulation of) discrete quantum walks on 
directed graphs is discussed. 

On the other hand, classically, multiple edges between 
a pair of vertices are equivalent to a single edge with a 
weight attached that modifies the probability of choos- 
ing that path to the next vertex. The generalization 
to weighted edges is straightforward, just replace the 
bits in the adjacency matrix by appropriately normalized 
weights. In the quantum case, different paths joining the 
same point can lead to interference effects, so multiple 
edges between two vertices are not equivalent to a single 
edge with a weight. We can also have complex weights 
for quantum walks, since they can be used to modify the 



amplitude of the wave function, rather than the proba- 
bility of choosing that edge. An alternative approach to 
multiple edges is to reduce it to the case of single edges 
by adding a vertex to each extra edge. This increases 
the size of the graph (iV) but keeps the degree of each 
vertex bounded (by N — 1). However, this reduction is 
not equivalent for all properties of interest on a graph, 
in particular for dynamical or flow properties, where the 
extra vertices and edges will slow down the progress of 
the random walk. We will not treat graphs with weights 
or multiple edges in this paper, though we note that this 
leaves interesting open questions about which types of 
graphs and properties are most amenable to analysis us- 
ing quantum versions of random walks. 

A. Discrete time classical walk 

A classical random walk on a graph G(V, E) can be 
performed as follows. Describe the location of the walk 
by an index x with < x < (N — I) running over the ver- 
tices in the graph. At each time step, a random number 
y between < y < (N — I) is generated, and the value 
of A xy examined to find out if edge e xy exists connecting 
the current location v x to v y . If the vertices are con- 
nected, the location index is updated to be y, (possibly 
only with some probability a, to tune the properties of 
the walk), and if not, the walk remains where it is at x. 
This procedure is repeated for as many steps as required. 
Note that we haven't bothered to exclude the case y = x. 
For large N , it isn't worth the extra computational time 
to check this condition for every new y generated. 

The probability P(x,t) of finding the walk at location x 
at time t is given by the difference equation representing 
one iteration of the walk, 

P(x,t+1) = P{x,t) + 

^Y,{A xy P(y,t)-A yx P(x,t)}.(l) 
y 

Since J2 y A yx — d x is the number of edges connected 
to vertex v x , there are various useful ways to rewrite 
this equation, especially if the values of d x are known in 
advance. 



B. Continuous time classical walk 

A continuous time classical random walk on this graph 
can be implemented in a similar way. We now have a 
hopping rate 7 per unit time to connected vertices, and 
the probability distribution can be written 

j f P(x, t) = 7 {A xy P(y, t) - A yx P(x, t)} . (2) 
y 

As written, this is clearly the continuous time limit of the 
discrete time version given in Eq. Q , with time rescaled 
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to t/N , and 7 = a. This formulation of a random walk is 
more useful where we have some type of analytic solution 
to Eq. |J2J) , such that the implementation has advantages 
over the discrete time random walk, or where it more 
closely matches a physical process we are trying to model. 

C. Resources required for implementation 

To quantify the resources required for implementing 
a random walk as an algorithm on a classical computer, 
observe that we can code the location index x in a binary 
string of size n = [log N~\ bits, exponentially smaller than 
the description of the graph (A) . We can assign the task 
of storing A to an oracle, and obtain the value of A xy 
by querying the oracle with the pair of numbers (x, y) at 
each step. This means the resources for performing the 
algorithm can be exponentially smaller than the descrip- 
tion of the graph. By assigning the description of the 
graph to an oracle, we are separating out the resources 
required to describe the problem from the resources re- 
quired to perform the algorithm. Specifying a problem 
or an answer to a problem can be highly non-trivial if it 
takes exponential resources to describe. A good exam- 
ple is solutions to Pell's equation, for which a quantum 
algorithm has been given by Hallgren 0] . 

This distinction is not merely formal. If we want to 
sample the properties of a random graph, we do not nec- 
essarily need to generate the whole description of the 
graph, we can just generate values of A xy on the fly as 
we sample them, allowing us to study the properties of 
structures larger than our computational resources can 
store. Also note that if we want to keep the algorithm 
small, i.e., 0(poly(n)), we cannot afford to increase the 
efficiency by remembering the values of (x, y) which we 
have already queried, since there are exponentially many, 
0(N ~ 2") of them even for one vertex. It follows that it 
can take an exponential number of steps to make useful 
progress on the graph, since the probability of making a 
move is ad x /N. If d x is 0(1), the chance of making a 
move is exponentially small. 

In the case of the continuous time random walk, the 
resource requirements are similar except we also have to 
account for the resources required to implement a suf- 
ficiently accurate numerical approximation of the con- 
tinuous time dynamics, possibly over exponentially long 
times, t ~ 0(N), or t ~ 0(N 2 ). 

D. More information about the graph 

The general graph we have chosen to define can thus 
be analyzed using a random walk algorithm with small 
requirements for memory, 0(poly(n)), but large require- 
ments for time, t ~ 0(poly(iV)). For many algorithmic 
and computational complexity purposes it makes sense 
to identify a class of problems for which the time re- 
quirements can also be small. This is the case for all 



the examples in the literature on which quantum walk 
dynamics have been analyzed. It means we either need 
to restrict ourselves to problems where N is not "expo- 
nentially large" , or else we need to know more about the 
graph, so that we can guarantee to move to a new lo- 
cation in 0(poly(n)) steps, rather than O(N). On the 
general graph we have no way of knowing ahead of time 
whether any pair of vertices is connected by an edge or 
not. But if we are given, say, a regular two-dimensional 
lattice with some of the edges missing, we can draw the 
random location index y from a much more restricted 
set (nearest neighbors, assuming they are labeled in a 
predictable way) and have a much higher success rate in 
finding the next move for the random walk. If we are do- 
ing a simple random walk on a line, cycle or hypercube, 
then we know exactly which edges are connected to any 
vertex, and the random choice is simply to pick which 
one to take. 

This reduction happens in general because we can en- 
code the description of a more regular graph in a more 
efficient way than the adjacency matrix. If the descrip- 
tion is exponentially smaller than the number of vertices 
N, then the description can be incorporated into the al- 
gorithm at no significant extra cost, and the oracle has 
no role to play in the analysis. Most of the quantum 
walks in the literature to date are on sufficiently regular 
structures that fall into this category. The main excep- 
tion is the "glued trees" graph of Ref . , which has half 
the edges in a random join down the middle of the graph. 
This still allows a quadratic reduction in the description 
of the graph, since the degree of the graph is fixed and 
0(1), so there are only O(N) edges to describe. However, 
the oracle is necessary here if we wish to use resources 
only 0(log N) for the algorithm itself. 

Though the description of the random walk process re- 
mains qualitatively the same in words, if the graph is now 
specified in some other, more efficient way than the ad- 
jacency matrix, Eqs. Q arid J2J) are no longer valid and 
must be replaced by versions using the appropriate graph 
specification reflecting the higher probability of making 
a move. It may be sufficient to adjust the value of a or 7, 
but often the equation for the dynamics of the random 
walk ends up looking rather different, once A xy has been 
replaced by a formula for generating the neighboring ver- 
tex indices, for example. 



III. QUANTUM WALKS ON GENERAL 
GRAPHS 

A quantum version of a random walk process must 
at least have the classical random walk as its classical 
limit. But this does not uniquely specify the approach to 
take. For a quantum evolution we require the dynamics 
to be unitary, and hence reversible, so the random ele- 
ment of the algorithm has to be catered for in other ways, 
for example, in the random result from a measurement 
with several possible outcomes. This gives us one route 
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to reproduce a classical random walk from the quantum 
dynamics, by applying a measurement at every step of 
the walk. Here we will first consider discrete time ver- 
sions of quantum walks, even though historically for a 
general graph such as we are considering, the continu- 
ous time quantum walk was presented first, by Far hi and 
Gutmann There are some subtleties to do with the 
continuous time case that will be best appreciated after 
the discrete time version has been described. 



A. Discrete quantum walks on general graphs 

Let us first consider the most naive way in which the 
classical discrete random walk on a general graph could 
be quantized. The location of the walk is represented by 
a state in a Hilbert space of size N, and may thus be in 
a superposition of positions. Intuitively, what we would 
like to do at each step of the walk is to distribute the wave 
function equally along all the connected edges from each 
vertex, instead of just picking one at random. However, 
the adjacency matrix is not unitary in general, so we 
cannot simply apply it as an operator to the particle state 
in a quantum version of Eq. QJ. We could try to make 
it unitary by introducing phases so the non-zero entries 
are e 1 ^ 1 " rather than unity, giving a quantum "adjacency 
matrix" A(Q) and 

\x,t+l) = | M ) + |:£{4Q)| y) i)-A«)M>} 
v 

= £{^?> + (i-f>4M>(3) 

y \. \ / J 

for the evolution. This only produces a unitar y o pera- 
tor under very restricted conditions, see Meyer who 
shows that this does not work for quantum walks on Eu- 
clidean lattices, and Severini [lJ,[Tj| for the more general 
case of graphs with directed edges. 

Instead, we need to model our quantum walk more 
closely on the algorithm for a classical random walk. This 
gives us a unitary process for any undirected graph. In 
the classical discrete random walk, we query the oracle 
with the pair of vertices (x, y) to obtain the value of A xy . 
For A xy = 1, the value of x is replaced with y and a 
new y generated for the next step. A quantum version of 
this procedure, with a superposition of positions instead 
of a single \x, y), needs to have the oracle receive and 
return the whole pair \x,y) along with a further qubit to 
hold the value of A xy . This is to ensure the operation is 
unitary, and also because we can't make and keep a copy 
of the state while passing \x, y) to the oracle. We can 
then perform a conditional swap operation depending on 
the value of A xy , to update the position of the walk. The 
quantum oracle O is defined by 

O\x,y,0) = \x,y,A xy ), (4) 



and the conditional swap S is defined by 

S\ x >y> A *y)-\ \x,y,A xy ) A xy = 0- W 

We aren't finished yet. We need to reset the value of the 
qubit holding A xy to zero, and generate a new value of 
y ready for the next oracle query. The first task can be 
done with a second call to the oracle, given the oracle is 
unitary and symmetric, and therefore its own inverse, 

G\x,y,A xy ) = \x,y,Q), (6) 

and since for the graphs we are considering, A xy = A yx 
it doesn't matter that some of the x, y pairs have been 
swapped. The value of y can be updated using any rea- 
sonable choice of unitary operation as a quantum "coin 
toss" C, that may be conditioned on the value of x if we 
want. Actually, the usual order is to toss the coin first, 
then apply the oracle-swap operations, so one step of the 
quantum walk is written 

m+i)) = osoc\m), (7) 

and t steps can be written 

\m) = (osocfim), (8) 

where |^(0)) can be the all zero state for walks starting 
at vertex zero. The only factor missing compared to the 
classical algorithm is a, determining whether we actually 
take the step (when allowed) or not. The quantum alter- 
native is to tune the properties of the walk by adjusting 
the coin toss operator C, if we want to keep the whole 
process unitary. This gives us more "knobs to twiddle" 
than we defined in the classical case, but in practice there 
are also many more ways to tune a classical random walk 
process than we have mentioned here. We could, for ex- 
ample, have a be vertex dependent, a(x), or even edge 
dependent, a(x, y), which is equivalent to weighted edges. 
Similarly, in the continuous time walk, 7 can be vertex 
or edge dependent. 

B. Resources required for discrete walk 

It is instructive to compare and contrast the resources 
needed for this quantum walk with the classical case. 
We have the same number of qubits as classical bits, i.e., 
2n + 1 , with n = [log N~\ . For the coin toss C on the y 
register, we should restrict ourselves to unitary operators 
that can be applied efficiently using 0(poly(n)) gates, or 
else we will have an exponentially inefficient algorithm. 
The simplest example would be to apply a Hadamard to 
each qubit in \y) individually, but the optimal choice will 
depend on the graph and problem to be solved 01 ■ The 
controlled swap S can be done efficiently as a sequence 
of Fredkin gates applied to each corresponding pair of 
qubits. 
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We need a quantum oracle rather than a classical or- 
acle, which is capable of processing and returning quan- 
tum superposition states. Given this quantum oracle, 
we need just two calls to it per step, whereas the clas- 
sical algorithm has to query the classical oracle N/d x 
times on average for each actual move made. This cer- 
tainly makes the quantum algorithm look more powerful 
for d x < N/2, since it can move up to m/2 steps for 
m oracle calls, compared to an average of md x /N steps 
for m oracle calls for the classical walk. But the quan- 
tum oracle is not the same as the classical oracle, which 
only returns one classical bit per query, so we need to be 
careful before drawing conclusions. 

If we insisted that the quantum algorithm must use 
the classical oracle, then it will need to retrieve the whole 
adjacency matrix in order to perform unitary operations 
that correspond to each possible pair of values of x and 
y that might be in superposition in the quantum state 
\ip(t)) of the walk. On the other hand, the classical walk 
could use the quantum oracle with the help of a qubit 
state preparation device, by querying it with a single pair 
of values \x, y, 0) and receiving \x, y, A xy ) in return, which 
could then be measured to obtain A xy with certainty. 

To illustrate further that the oracle-algorithm separa- 
tion is fairly arbitrary in this situation, we can redefine 
the oracle to be the whole conditional swap operation 
OSO. Classically this just means we give the pair (x,y) 
to the oracle, and receive in return either (x, y) if there 
is no edge connecting those two vertices, or (y, x) if there 
is a connecting edge. Since the numbers are classical 
bits, we could inspect their values and detect whether 
the swap took place, but we don't need to know, we just 
toss a new value for y and continue; we no longer need 
the single bit for the value of A xy . The quantum version 
of this oracle is 

oso= J2 \v> x )( x >v\ + E \ x >vX x >v\> (9) 

where we also need one less qubit, since the value of A xy 
is not returned. 

On a more basic level than resource counting, we have 
not shown that this definition of a discrete quantum walk 
does anything usefully different from a classical random 
walk in the general case. However, the specific cases of 
coined quantum walks in the literature are clearly deriv- 
able from it by encoding the simpler graph more effi- 
ciently For example, the walk on a hypercube is obtained 
by noting that there are exactly n choices of connected 
vertices rather than a possible (N — 1), and shrinking the 
size of the \y) qubit register accordingly to [logn]. The 
oracle operation is redundant (unless the labeling of the 
vertices has been randomized) since the value of y is used 
to specify which bit of the value of x to update to obtain 
the new vertex index. Useful classical random walk algo- 
rithms are also usually highly tuned to suit the problem, 
so we shouldn't necessarily expect that a proof of quan- 
tum superiority over classical would be forthcoming in 
the general it will depend on the particular choice 



of coin toss operation |18| and the particular property of 
the graph we are trying to determine. It has already been 
shown that mixing times on a hypercube are not faster 
than classical 0], whereas the hitting times are, at least 
for the equivalent classical random walk process |2jj ■ 



C. Continuous time quantum walks on general 
graphs 

The continuous time quantum walk on a general graph, 
as presented by Farhi and Gutmann 0, simply uses 
the adjacency matrix A, which is Hermitian for an undi- 
rected graph, to form the Hamiltonian for the evolution 
of the quantum state, 

4<*w*)> = £<*i h mw(*)>> ( 10 ) 

y 

with H = 7 A, and solution 

m)) = e^ At \m)- (ii) 

Comparing this with Eq. J2J shows that the second 
term, necessary for conservation of probability, is miss- 
ing. Since we only need a Hermitian operator here, and 
since A xy = A yx guarantees this, we are free to examine 
this quantum evolution as well as that obtained in more 
direct analogy to Eq. (J2J, 

y 

- (y\A\x)(x\m)}- (12) 

For graphs where all the vertices are of the same de- 
gree d, the second term only introduces an irrelevant 
global phase, therefore it makes no difference to observ- 
able quantities if it is omitted |2l| . The Hamiltonian in 
Eq. I|12l) becomes H = 7(A — rfl), and the solution to 
this can be written 

|^(i)> = e -^ A ~ dl )*|</;(0)). (13) 

Since A commutes with the identity, the two terms in 
the exponential can be factored, giving 

|^(t)) = e -* 1 ' A *e < T , * 1 *[V(0)>, (14) 

which is the same as Eq. Ill 1^ apart from the phase e* 7d -"-*. 
For graphs of general degree, however, the dynamics with 
the second term included will have a different evolution 
|22| . So far, only graphs of fixed degree have been studied 
in any detail in the literature so the difference between 
the two versions has not been explored. 

Equation Hll|) looks remarkably similar to the discrete 
case, Eq. (|SJ|, both can be written in the form 

\m) = (w)*wo)>, as) 
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with unitary operator U — OSOC for the discrete time 
walk and U = e~' nA for the continuous time walk. But, 
unlike the classical case, where Eq. (J2J is the limit of 
Eq. Q as the time step goes to zero, in the quantum case 
the similarity is deceptive, the discrete and continuous 
walks have Hilbert spaces of different sizes, Hn in the 
continuous case and Hn ® Hn <E> Hi in the discrete case. 
Moreover, though the H2 component returns to |0) at 
the end of each step and thus factors out, the discrete 
walk evolution does not remain in an Hn subspace (apart 
from trivial special cases), so this is not a simple route 
to showing how the two versions are related. 

D. Resources required for continuous walk 

For algorithmic purposes we need to encode the walk 
into n-qubit registers, and evaluate Eq. (II 1|) for our cho- 
sen values of t and |-0(O)). Childs et al. 1], in their quan- 
tum walk algorithm, showed how to accomplish this for 
a graph of degree 0(1), where their oracle held a more 
efficient description of the graph than the adjacency ma- 
trix. Here we will describe their method as applied to a 
general graph, using the same oracle O as in the discrete 
time quantum walk defined by Eq. Q. 

We take a qubit register of size n to hold the vertex 
index indicating the current location of the walk, and 
two ancillary registers, one also of size n qubits and the 
other a single qubit. The Hilbert space is thus of size 
Hn ® Hn <8> Hi , the same size as for the discrete time 
quantum walk in Sec. IIII Al First we will describe what 
goes wrong if we try to avoid introducing a "consistent 
coloring" of the edges. Define the set of unitary operators 
V y that write the value of y into the second register, 

V y \x,z,b) = \x,z®y,b), (16) 

where © specifies bitwise addition, and we will take the 
initial values of z and b to be zero. Note that V y is its 
own inverse. After applying V y , we can apply the oracle 
O to place the value of A xy in the last qubit. Since the 
actual input can be a superposition of any of the N vertex 
indices x, the operation OV y can produce a superposition 
of all possible A xy for a given y. Next we apply the 
Hermitian evolution that does a conditional swap on the 
first two registers, and keeps only the terms in the graph 
with A xy — 1, as indicated by the final qubit, 

T\x,y,l) = \y,x,l) 

T\x,y,Q) = 0. (17) 

This can also conveniently be written 

T\x,y,A xy ) = A xy \y,x,A xy ). (18) 

Childs et al. pj show that e _lTt can be evaluated effi- 
ciently using (9(poly(n)) gates, and give a circuit diagram 
for this. We will not repeat these steps here. Next we 
must reset the extra registers to zero, in order to keep 



the evolution of the walk in the subspace of the first n 
qubits. Removing the value of A xy from the last qubit is 
accomplished with another application of the oracle O, 
but, because the first two registers have been swapped, 
we now need to apply V x to get what we want 

V x |»,x,0) = |y,0,0>. (19) 

However, since V x 7^ V y , this won't give us a Hermitian 
operator overall. This is where the consistent coloring 
comes in. The edges are given labels (colors) such that 
all vertices have at most one edge with each color. Then 
the operators are labeled by color instead of by the vertex 
they connect to, and the quantity that is added to the 
second register is y c (x), the index of the vertex connected 
to v x by color c, 

V c \x,z,b) = \x,z®y c (x),b) (20) 

Again, V c is its own inverse, but now it also provides 
the required operation when applied after T to reset the 
auxiliary registers 

Vc|y, x, 0) = \y,x® y c (y),0) = \y, 0, 0). (21) 

So 

H = 7 VcOTOV c (22) 

is a Hermitian operator that enacts the dynamics in 
Eq. ((TTJl . This can easily be checked by calculating 
(y, 0, 0|H|a;, 0, 0} and showing that it equals jA xy . Also, 
since H|x,0, 0} = J2 c A x, yc (x)\y c (x),0,0), provided the 
walk starts in the subspace spanned by the first set of n 
qubits, it remains in this subspace at the end of the walk. 

A consistent coloring is always possible using one more 
color than the maximum degree of the graph 23]. Here, 
since we only know the maximum degree of the graph 
is bounded by (N ~ 1), we have to assume a coloring 
using N colors. This at least makes the coloring trivial to 
implement, since a coloring that works for the complete 
graph (all vertices joined to all others) will work for any 
other graph. Such a coloring can be generated by using 
c = (x + y) (mod N) to create the color label c for edge 

&x y • 

To simulate a Hamiltonian given as a sum of terms as 
in Eq. (|22|l . we make use of the Lie product formula, 

e -i{H 1 +H 2 + ---+H N )t 

= lim f e -i^t/j e -iH 2 t/j e -iH„t/j\ 3 t (23) 

only we use a finite value of j chosen to give the desired 
accuracy. Lloyd [24j showed this is possible using j ~ 
0(poly(A^)), where N is the number of terms in the sum. 
We will thus have a total of Nj ~ 0(poly(iV)) operators 
of the form V c Oe _i7Tt / J OV ' c to perform. 

For the specific case treated by Childs et al, this 
scheme provides an efficient algorithm in the sense that 
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it requires only 0(poly(log N)) gates to implement be- 
cause they know the degree of the graph is O(l) [in fact 
d = 3 apart from two special vertices with d = 2], so 
the number of terms in the sum over colors in the Hamil- 
tonian only needs to be 0(1). However, in the general 
case we are considering, where the degree of the graph is 
bounded only by N — 1, we have N terms in our sum over 
colors, so we end up with 0(poly(A)) gates rather than 
0(poly(log A)). This should be seen as a consequence of 
the general nature of the graph we defined, see Sec. Ill Ul 
rather than a failure to find an efficient algorithm. In 
the discrete time quantum walk, we may expect to need 
0(poly(iV)) time steps, and thus 0(poly(iV)) quantum 
gates too. 

The correspondence between the quantum computer 
programs for the discrete and continuous time walks goes 
deeper than this. In order to provide a finite precision 
implementation of the continuous time walk on a quan- 
tum computer, we have broken down the continuous time 
Hamiltonian evolution into discrete steps thus, 

N-l 

\m) = II (VcOe'^ Tt ^OV c y |V(0)>. (24) 

If we insert 1 = W-iVjv-i at the start and rewrite 
slightly as 

N-l 

\ip(t))=V N -i I] (Oe-^'WVcVc-^ Viv-ilV(O)), 

c=0 

(25) 

where we take c — 1 (mod N) so when c = we have 
c — 1 = N — 1, then the correspondence with the discrete 
walk, \i/)(t)) = {OSOCf\ip{{))) , Eq. © becomes clear. A 
pair of V c operators take the place of the coin toss C, and 
the evolution e~ nTt ^ replaces the controlled swap S. 
These operators even perform the same function as their 
discrete counterparts: choosing a new value of y, and 
swapping the x and y values when A xy — 1, respectively. 
The pair V c V c -i varies from step to step as c varies, so 
this evolution is using a variable coin toss operator as well 
as a different sort of swap operation, but the structure is 
essentially the same. 

As already noted, there are Nj steps to this discretiza- 
tion of the continuous time quantum walk. Each unitary 
operator besides the oracle can be performed efficiently, 
using (9(poly(log N)) gates, so each step uses roughly the 
same number of qubits and gates as a single step of the 
discrete time walk. Moreover, the discrete and continu- 
ous versions of the quantum walk are likely to be run for 
roughly the same number of steps. There is thus no es- 
sential difference between the algorithms in terms of the 



number of oracle calls or the number of quantum gates. 

IV. SUMMARY AND DISCUSSION 

Quantum walk algorithms give us a new type of algo- 
rithm to run on our quantum computers (when they are 
built), but physical quantum walk processes don't give us 
a new type of quantum computer architecture unless the 
walk takes place in an efficiently encoded Hilbert space 
rather than a physical position space. 

The discrete and the continuous time quantum walk 
algorithms on a general graph have similar structures 
when implemented on a quantum computer, with sim- 
ilar efficiency and resource requirements. This does not 
prove their equivalence in any algorithmic sense, but it 
does provide supportive evidence for the two formulations 
having the same computational power. 

All the quantum walk algorithms discussed here de- 
pend on the graph being undirected, because they rely 
on the adjacency matrix being symmetric. So there are 
parallels with algorithms based on the quantum Fourier 
transform, which can solve the hidden subgroup prob- 
lem for Abelian groups and some extensions, but not in 
general for non- Abelian groups. 

The question of how the discrete time quantum walk 
relates to the continuous time quantum walk remains 
a key open question for the theoretical development of 
quantum walk algorithms. A more practical challenge 
is to look for quantum measurements that extract the 
kind of information typically gathered from a classical 
random walk process. These are often average or aggre- 
gate quantities that may conceivably be accessible with a 
single measurement, even though direct measurement of 
the quantum state of the quantum algorithm would yield 
only one value for the position of the walking particle. 
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